Abstract. We have mapped the dark cloud LI 83 in the far- 
infrared at 100/^m and 200//m with the ISOPHOT photome- 
ter aboard the ISO satelHte. The observations make it possi- 
ble for the first time to study the distribution and properties of 
the large dust grains in LI 83 without confusion from smaller 
grains. The observations show clear colour temperature vari- 
ations which are likely to be caused by changes in the emis- 
sion properties of the dust particles. In the cloud core the far- 
infrared colour temperature drops below 12K. The data allow 
a new determination of the cloud mass and the mass distribu- 
tion based on dust emission. The estimated mass within a ra- 
dius of 10' from the cloud centre is 25 Mq. We have mapped 
the cloud in several molecular lines including DC0h-(2-1) and 
H^'^COh-(I-O). These species are believed to be tracers of cold 
and dense molecular material and we detect a strong anticorre- 
lation between the DCOh- emission and the dust colour temper- 
atures. In particular, the DCOh-(2-1) emission is not detected 
towards the maximum of the 100/J,m emission where the colour 
temperature rises above 15 K. The H^^COh- emission follows 
closely the DCOh- distribution but CO isotopes show strong 
emission even towards the lOO/im peak. Detailed comparison 
of the DCOh- and C^^O maps shows sharp variations in the rela- 
tive intensities of the species. Morphologically the 200/im dust 
emission traces the distribution of dense molecular material as 
seen e.g. in C^^O lines. A comparison with dust column density 
shows, however, that C^^O is depleted by a factor of ~L5 in 
the cloud core. We present results of R- and B-band starcounts. 
The extinction is much better correlated with the 200 /im than 
with the 100/zm emission. Based on the 200 /im correlation at 
low extinction values we deduce a value of '^17'" for the visual 
extinction towards the cloud centre where no background stars 
are observed anymore. 
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1. Introduction 

The distribution of material in dense interstellar clouds can 
be traced by either molecular line emission or the continuum 
emission from dust particles. The analysis of line spectra is 
complicated by radiative transfer effects. The emission of ra- 
diation depends on local properties like density and kinetic 
temperature. Based on these parameters, the observed line in- 
tensities can be predicted with radiative transfer modelling. In 
practice, the density structure of the sources is, however, un- 
known and e.g. small scale density and velocity variations can 
seriously affect the observed intensities. Secondly, the ratio be- 
tween the abundance of the tracer and the total gas mass is 
only known approximately and it can change even within one 
source. The chemical reactions depend on density and tempera- 
ture (Millar et al. |1997l Lee et al. |1996|) and therefore the chem- 



istry of a dark cloud cores will differ from the chemistry in the 
surrounding warmer material. 

Molecular abundances depend on the strength of the ultra- 
violet radiation field via processes such as photodissociation. 
The radiation is weakened by dust extinction and by the self- 
shielding of the molecules and this introduces radial abundance 
gradients. Selective photo-dissociation can lead to strong abun- 
dance variations between different CO isotopes and can even 
affect the excitation of the molecules (Warin et al. 1996l The 



abundances change even as a function of time (Leung et al. 
1984; Herbst & Leung 1989, 1990; Bettens et al. 1995) and 
the observed abundances have been used as a 'chemical clock' 
to determine the evolutionary stage of sources (e.g. Stahler 



1984 Lee et al. |1996| ). Observations of e.g. TMC-1 show clear 
:arbon bearing and oxygen bearing 
1997). This can be compared with 



anti-correlation between carbon bearing and oxygen bearing 
molecules (Pratap et al 



and long carbon chain molecules are produced preferentially 
at early times, while SO abundance rises only much later. 

The role of dust grains in chemical evolution is important 
but not yet well understood. In dark cores some molecules can 
freeze onto the dust grains and this causes the under-abundance 
of the molecules in the gas phase. CO abundances can decrease 
by more than a factor of ten (Gibb & Little |1998| ). In the cold 
cores, the abundance of deuterated species, e.g. DCOh- are en- 
hanced. This fact has been used to identify condensations in the 
earliest stages of star formation, before the birth of protostars 
again increases the temperatures (e.g. Loren et al. 1990). The 



abundance ratio [DCOh-]/[HCOh-] depends also on the abun- 
dance of free electrons and can therefore also be used to deter- 
mine this quantity (e.g. Guelin et al 
Wootten 



1977 



1982; Anderson et al. 1999l. The electron abundance 



1982 



Watson 



1978 



is an important parameter since, through ambipolar diffusion, 
it affects the timescales of physical cloud evolut ion an d even- 
tually the rate of the star formation (Caselli et al. 1998 ). 

The analysis of far-infrared dust continuum observations 
is in principle more straightforward. The emission is opti- 
cally thin and the intensity directly related to the dust column 
density. At wavelengths 100-200/im the radiation is emitted 
mainly by the classical large grains, with a possible small con- 
tribution at 100 fim by the 'very small grains'. The large grains 
heated by the normal interstellar radiation field settle to a tem- 
perature of ~17 K and the maximum of the spectral energy dis- 
tribution is close to 200/im. The temperature of the small grains 
fluctuates and most of the radiation is at sh orter w avelengths 
(A < lOO^im). In the model of Desert et al. ( |199C| ) some 10% 
of the infrared emission at 100/im comes from very small dust 
grains. 

The large grains are well mixed with the gas phase (Bohlin 



models (e.g. Bergin et al. |1995[ |1996|) which show that CS gt al. and the FIR emission can be used to trace the gas 
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distribution. The dust emission corresponding to a given col- 
umn density varies depending on the composition and tem- 
perature of the dust. The far-infrared colour temperature has 
been observed t o drop towards the centre of many dark cores 
(e.g. Clark et al. 1991 ; Laureijs et al. 1991 ) and even cirrus type 



clouds (Bernard et al. 1999). It is clear that the variations can- 
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not be entirely due to a reduced radiation field, i.e. the ob- 
served colour temperature variations indicate a change in the 
grain properties. In the cloud cores, the dust particles may grow 
through coagulation, and in this case the associated change in 
the grain emissivity could explain part of the colour tempera- 



ture variations (e.g. Meny et al. 200C; Cambresy et al. 2001) 



Changes in the dust properties cause uncertainties to the col- 
umn density estimates. These are, however, smaller than those 
associated with the analysis of optically thick molecular lines. 

Most of the listed problems can be addressed by comparing 
observations of molecular line emission with far-infrared dust 
emission. In this paper we present such a study of the dark 
cloud L183. 

1.1. L183 

L183 (L134N) is a well known dark cloud belonging to a high 
latitude cloud complex (Galactic coordinates {I, b)={6.0, 36.7)). 
The dista nce of th e clou d has been estimated to be lOOpc (Mat- 
tila 1979, Franco 1989) although a distance of 160pc has also 
been used (Snell 1981). We adopt 6?= 100 pc. Extinction towards 
the centre of L183 is at least Ay ~6 (Laureijs et al. 1991) and 
may exceed Ay —10. The cloud core is surrounded by an ex- 
tended low extinction (Ay <1) envelope. 
L183 was observed in CO(l-O) and ^^€0(1-0) by Snell 



( 1981 ), in NH3 by Ungerechts et al. (198C) and in sever al lines 
includi ng Ci^O(l-O) and H13cOh-(1-0) by Swade ( |l989a ; 
1989t). The ^"^CO line map shows a cometary shape with the 



head of the globule pointing north in equatorial coordinates 
(~galactic east). This is also the densest part of the cloud. The 
cloud has a sharp edge towards the north, but in the south- 
ern part the strong ^■^CO and especially the ^^CO emission 
continues in the tail that points towards L1780. Together, the 
region containing clouds L134, L169, L183, L1780 is some- 
times called the LI 34 complex. NH3 and C3H2 emission is 
concentrated into the head region and ammonia emission is re- 
stricted to a very narrow ridge running from south to north. 
The usual centre position given for the cloud, (RA=15'*51™30'* 
Dec=-2°43'31", 1950.0) is close to the centi-e of the ammo- 
nia emission where density reaches a few times 10^ cm 
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(Swade 1989b ). In recent years the cloud has been observed in 



a number of molecular lines (see e.g. Swade & Schloerb 1992; 



Turner et al. 1999, 2000 and Dickens et al. 2000 and references 
therein). 

The infrared emission from LI 83 was studied by Laureijs 
et al. ( |1991[ |1995| ) based on the IRAS data. The 60^m emission 
is detected only in a narrow layer surrounding the cloud core. 
Based on the upper limit of the I(60/im)/I(100/im) ratio, an up- 
per limit of 15 K was derived for the colour temperature assum- 
ing a ly^ emissivity law. However, at 60/im, most of the emis- 
sion comes from small dust particles, and the colour tempera- 
ture is not directly related to the physical grain temperature. In 
the cloud core the lack of 60/im emission was interpreted as ev- 



idence for mantle growth. Ward-Thompson et al. (1994; 200C) 



Lehtinen et al. (2000) have studied the properties of the sources 
using far-infrared observations made with ISO. 

In this paper we study the dark cloud L183 using 100 /im 
and 200 /im maps and new molecular line observations. The 
far-infrared observations were carried out with the ISOPHOT 
instrument aboard the ISO satellite. Both 100/im and 200/im 
are dominated by emission from classical large grains and this 
makes it possible, for the first time, to study the variations in the 
properties of the large grains without confusion from smaller 
size dust populations. The far-infrared observations enable us 
to determine the mass distribution of L183 independently from 
the molecular line data which may be affected by chemical 
fractionation and other systematic effects. The dust distribu- 
tion will be compared with molecular line observations which 
include new ^^CO, ^^CO and C^^O observations that cover the 
entire cloud LI 83. The correlations with dust column densi- 
ties will be used to determine the degree of C^^O depletion in 
the core of LI83. New DCOh-(2-1) and H13cOh-(1-0) map- 
pings have been made with higher spatial resolution than what 
were previously available. DCOh- is known to be a tracer of 
cold gas and we will determine the correlations between DCOh- 
and dust emission, especially the FIR colour temperatures. We 
will also study the variations in the intensity ratio of DCOh- 
and H^'^CO-H. Finally, we present results of R- and S-band star 
counts and a comparison between the extinction and the far- 
infrared emission. 



2. Observations and data processing 

2.1. FIR observations 

Far-infrared observation of L183 were done with the ISOPHOT 
photometer (Lemke et al. 1996) aboard the ISO satellite 
(Kessler et al. 1996). Two 30' x 30' maps were made in the 



have made sub-millimetre continuum observations of the cloud 
and have detected two point sources close to the cloud centre. 



PHT22 staring raster map mode with filters C_100 and C_200. 
The reference wavelengths of the filters are 100 /xm and 200 /im 
and the respective pixel sizes ~44" and 89". 

Additional narrow strips running from NW towards the 
cloud centre were observed at 80/im, 100/im, 120/im, 150/im 
and 200/im. The analysis of these observations will be pre- 
sented by Lehtinen et al. (in preparation). 

Comparison of the ISOPHOT 100/im observations with 
the 100 /im IRAS ISSA maps showed peculiar scatter with two 
separate branches in the ISOPHOT vs. IRAS relation. The cali- 
bration of the ISOPHOT observations was done with Fine Cali- 
bration Source (FCS) measurements that were done before and 
after the actual mapping and a linear interpolation was per- 
formed to determine the detector responsivity at each map po- 
sition. The detector signals did not, however, stabilize during 
FCS measurements and this leads to some uncertainty. A 5% 
change in the relative responsivities of the two FCS measure- 
ments was found to be enough to minimize the scatter rela- 
tive to the IRAS map. The average responsivity over the map 
was not changed. The resulting relation between ISOPHOT and 
IRAS ISSA surface brightness values is shown in Fig. |l]. Be- 
cause of the large temperature variations the relation between 
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Fig. 1. The relation between the 100/im ISOPHOT and IRAS 
ISSA surface brightness values. The ISOPHOT time-ordered 
data have been tilted 5% in order to obtain lowest dispersion 
when compared with IRAS values (see text). The slope of the 
fitted line is 1.09. No zodiacal light has been subtracted from 
the ISOPHOT values. 



100/im and 200/im surfaces brightness values shows a signifi- 
cantly larger scatter (see Fig. This is also decreased by the 
5% adjustment of the 100/im responsivities. 

The 100/im map was finally rescaled to the DIRBE sur- 
face brightness scale by multiplying the values with 0.85. 
The scaling was established by correlating DIRBE, IRAS and 
ISOPHOT surface brightness values in the area. Details of the 
procedure are given in Lehtinen et al. (in preparation). The 
200 fj,m observations were calibrated using the PCS measure- 
ments which were of good quality. The uncertainty in the ab- 
solute calibration of the 200 /im map is expected to be below 
20% (see Garcia-Lario 2000 ). The final far-infrared maps are 
shown in Fig. |[ 



2.2. Molecular line observations 

The molecular line observations were made with the SEST 
telescope during two session in February 1998 and February 
1999. During the first observing period, the cloud was mapped 
in the lines 1^00(2-1), 13CO(1-0) "C0(2-l), C34S(2-1) and 
C^*O(l-0). Due to a problem with the receiver, the C^'*0(2- 
1) map remained uncompleted. In February 1999 observations 
were made in the DCOh-(2-1) and H13cOh-(1-0) Hnes. 

During observations, the pointing was checked every two 
or three hours by observing bright SiO maser sources. The es- 
timated pointing uncertainty is less than 5". Most observations 
were made in the frequency switching mode with a frequency 
throw 6 MHz. The chopper-wheel method was used for calibra- 
tion and observed intensities are given in T^ units. The typical 
rms noise levels are given in Table |l[ 
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Fig. 2. Observed 100 /im surface brightness as the function of 
the 200 /im values. The diffuse background determined from 
regions with lowest surface brightness has been subtracted 
from the observed values and the 100/im data are convolved 
to the resolution of the 200/im observations.. The solid lines 
correspond to dust colour temperatures 12.5 K and 15.0 K. 
The dashed line shows the least squares line fitted to points 
5(200 /im)<30MJysr-i 

Table 1. Summary of observed molecular emission lines. The 
columns are: (1) transition, (2) beam size, (3) spectral resolu- 
tion (4) number of observed positions and (5) the average rms 
noise. 
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FWHM 
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< ctrms 
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(K) 


'^CO(2-l) 


23 


0.054 


83 


0.17 


^='CO(1-0) 


46 


0.114 


83 


0.069 


^=^00(2-1) 


24 


0.057 


532 


0.14 


C^**O(l-0) 


46 


0.114 


533 


0.10 


H^3CO+(1-0) 


57 


0.144 


174 


0.040 


DC0-H(2-1) 


35 


0.087 


176 


0.061 



The ^^CO, ^'^CO and C^^O maps cover practically the en- 
tire cloud LI 83 as seen in optical extinction. In the southern 
part, a fairly strong molecular line emission continues up to 
the edge of the map. The H13cOh-(1-0) and DCOh-(2-1) maps 
are smaller in size, but particularly in the case of DCO-n, they 
cover the whole region of significant emission. At the map 
boundaries, the DCOh- intensity has dropped from a peak value 
Ta wl.lKto -0.25 K or below. 

2.3. Star counts 

The optical extinction in the LI 34 cloud complex was deter- 
mined by means of star counts using a blue (IlaO + GG385 
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Fig. 3. a-b) The 100/im and 200/im I SOPH OT maps of the cloud L183. The cross indicates the position of the continuum source 
detected by Ward-Thompson et al. ( |l994| The (0,0) position is RA=15''51™30^ DEC=-2°43'31". In the case of 100/im the 
contours are at 15 and 20MJysr^^ and for 200/im the contour levels are 30, 50 and 70 MJysr^^. The small squares at the 
lower left corner indicate the pixel sizes, 44" for 100 fxm and 89" for the 200 fim map. c) The map of dust colour temperature 
calculated from the 100/im and 200/im ISOPHOT observations. The contours are at intervals of 1 K starting with 12 K. The 
effective resolution is 100". d) blue image from the Digitized Sky Survey. 



filter) and a red (IIIaF + RG630) ESO Schmidt plate (Plate 
Nos.3713 and 3714) which were obtained on 16 April 1980. 
Star counts were performed using a reseau size of 2.8' x2. 8'. 
The total number of stars down to th e limiting magnitude of 
each plate was counted (see Bok 195f for the method). 



In order to derive the extinction, the function N{m) is 
needed, i.e. the number of stars per sq. degree with magnitude 
< m in a transparent comparison area close to the dark nebula. 
The tables of va n Rhij n ( |1929 ) have often been used for this 
purpose (see Bok 1956 ). In the present case it was better to use 
the RGU photometric catalogue of Becker and Fenkart ( 1976| ) 
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Fig. 4. The extinction maps of the LI 34 dark cloud complex 
derived from starcounts in the R- (lower frame) and B-bands 
(upper frame). The area mapped with ISOPHOT is indicated 
with a square. Individual clouds LI 34, LI 83 and LI 69 are in- 
dicated in the figure. The contour levels are at 1, 2 and 4™. 



for Selected Area 107 which is located close (1=5.7 deg, b= 
4L3 deg) to the LI 34 complex. The R and G band star numbers 
for SA 107 were extracted from the compilation by Bahcall et 



al. (1985). These N{m) curves were extrapolated beyond the 
limiting magnitudes R=\l .5 and 0=18.5, and were corrected 
to account for the somewhat different Galactic latitude (6 = 36 
deg) of the LI 34 compl ex by using the tabulated results of Bah- 
call and Soneira ( 1980 ) for their Galaxy model in the B and V 
bands. 

An essential parameter of the logN{m) vs. m caUbration 
curve is its slope at the limiting magnitude of the star counts. 
In our case this slope turned out to be 0.22 at Gum = 20. 1 mag 
and i?ii„i=19.3 mag, respectively. We note that the G and R 
passbands of the Basle RGU photometric system with effective 
wavelengths of 463 and 638 nm are close enough to the wave- 
length ranges of our blue (380^80 nm) and red plate (630- 
690 nm). 

The resulting extinction maps are shown in Fig. Q 



3. Observational results 

3.1. FIR maps 

The main structure seen in the far-infrared maps (Fig. ||a-b) is 
a wedge pointing northwards. The intensity is highest close to 
the head of the wedge, and drops rapidly towards the north and 
more slowly towards the south. There are clear differences in 
the distributions at 100/im and 200/im. There are two broad 
200/im maxima, one close to the optical extinction peak and 
previously known mm- and sub-mm continuum source (Ward- 
Thompson et al. 1994) and another ^^5' towards SE. In the 



100/im map only the latter peak is visible. The relative strength 
of the 200/im emission is larger in the north, west and south- 
west while the lOO/xm emission is stronger in the southeast. 
In the west the 100/im emission seems to extend further out, 
although the surface brightness is already close to the back- 
ground level. 

The 200/im emission is caused entirely by the large clas- 
sical grains and at 100/im the contribution from small grains 
is expected to be small, ^10%. Dust temperature map was de- 
rived assuming that large grains emit according to a modified 
black body radiation v^B{v^ Td). The background correspond- 
ing to the lowest surface brightness values within the mapped 
area was first subtracted. Average 100/im surface brightness 
values were calculated corresponding to each 200/im measure- 
ment. Values were colour corrected and the dust temperature 
was obtained by the fitting of the modified Planck curve. Since 
the colour correction depends weakly on the assumed temper- 
ature the procedure was repeated once using the temperature 
derived from the initial fit. 

The dust temperature in Fig. |]: has a well defined minimum 
close to the position of the continuum source. The minimum 
temperature in the 200/im core is below ^ 12K while in the 
north and west the temperature rises rapidly to '^16K. There 
is no temperature maximum at the 100 /im emission maximum 
and the position does not differ in any way from its surround- 
ings. The 5% tilting of the 100/im map (Sect. 2. 1 ) changed the 
temperatures by up to ^1 K at the eastern and western borders 
(the first and last scan lines in the map). As the average respon- 
sivity was not changed the temperature in the centre, i.e. in the 
cold core, was not affected. 

The correlation between 100/im and 200/im surface 
brightness values (Fig. ^ shows considerable scatter. The scat- 
ter was somewhat reduced by the correction applied to the 
relative responsivity of the 100/im FCS measurements (see 



Sect. 2.1 ). A further modification of the responsivities would 



have very little effect on the scatter, which is caused by true 
variations of the dust properties. The variations are reflected in 
the dust temperature distribution. 

A linear least squares fit to the data points in Fig. ^ 
gives a relation 5fit(200/im)= (0.14±0.01) S(100/im) + 
(1.30±0.06). The fitting procedure takes into account the 
uncertainty in both variables. The statistical error of the 
intercept in particular is small compared with the uncer- 
tainty due to calibration and the background subtraction. 
Data at positions where the background subtracted values of 
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Fig. 5. Distribution of AS'(200 /^m). AS'(200 fim) is the differ- 
ence between the observed 200 fim surface brightness and the 
least squares prediction calculated from the 100 fim values. 




EO 10 -10 



ia [arcmlu] 

Fig. 6. 200 /im surface brightness (contours) and the calculated 
optical depths at the same wavelength. The contours start at 
15 MJy sr^^ and are drawn at intervals of 15 MJy sr^^. 



S{2QQ fim) were below SOMJysr ^ (i.e. outside the cloud 
core) give a relation Sm(200 ^im)= (0.22±0.01) S(100/im). 
The limit of 30MJysr~^ at 200 /im corresponds approxi- 
mately to the 50MJy sr^^ contour in Fig. ^ where no back- 
ground subtraction was made. According to Fig. ^ the ex- 
cluded area corresponds also to the region with highest ex- 
tinction, Ab ^4"\ Based on the derived dependence we cal- 
culate AS'(200/xm)=S'obs(200yLtm)-S'fit(200Aim), which is the 
difference of the observed 200 fim surface brightness and the 
predicted value based on the 100/im data. The resulting map 
of AS'(200 Aim) (Fig. |]) shows the 200 /im excess that can be 
due to either cold dust or dust with enhanced emissivity at the 
longer wavelength. 

The 200 /im optical depth was calculated from the observed 
intensities as t(200 /im)=5/B(r). Here B{T) is the black 
body intensity at temperature T as read from Fig. The re- 
sulting map, which shows the column density distribution of 
large dust grains, is rather similar to the 200 /im intensity dis- 
tribution. Due to its low colour temperature, however, the main 
core is much more pronounced. The region of the 100 /im peak 
is correspondingly weaker and is visible only as an extension 
of the main core. Fig. || shows the 200 /im surface brightness in 
relation to the 200 /im optical depth. While the SB core is also 
visible as a separate intensity peak in the 200 /im map, the large 
grains are clearly concentrated in the NW core. A similar rela- 
tion is seen between the 200 /im emission and the dust colour 
temperature. In the temperature map the SB peak is completely 
invisible. 

The cloud mass was calculated from the 200 /im optical 
depth using the value cth = 2.5 • lO^^'^cm^ for the effective 
dust cross section per hydrogen atom at 200 /im (Lehtinen et 



the value of cth and use of the formula given by Boulanger 
1996|, (TH = 10"25(;^/250/im)-2cm-2, would result in 



et al. 

some 50% higher mass estimates. We note that the latter an 
value is for diffuse medium, while the Lehtinen et al. value is 
for a dense globule similar to LI 83. 

As shown in Fig. ||, there is a clear 200/im excess in the 
centre of L183. This could be due to a decrease in the physical 
temperature of the dust grains induced by the decreasing inten- 
sity of the radiation field, or by increased overall FIR emissivity 
of the grains. Another possibility would be a change in grain 
properties, i.e. increased emissivity at the longer wavelength, 
leading to a change of the dust emissivity index. Both explana- 
tions are related to the increasing density in th e cloud centre. 
Both theoretical studies (e.g. Ossenkopf |1993| ; Wright [l987| ) 
and rece nt PRO NAOS observations (Bernard et al. 1999; Step- 
nik et al. 2001 ) indicate that dust emission properties are likely 
to change in dense and cold clouds. This could be the result of 
grain growth, or possibly the formation of grain aggregates. In 
a forthcoming paper (Juvela et al., in preparation) we will study 
the far-infrared dust emission using radiative transfer models. 
Further discussion about the dust temperatures and changes in 
the dust emissivity law will be deferred to that paper. 

3.2. Distribution of molecular line emission 

Figure ^ shows the distributions of the integrated antenna 
temperature in lines i2cO(2-l), i3co(2-l), Ci**O(l-0), 
H13cOh-(1-0) and DC0h-(2-1). 

The CO lines follow the general morphology seen in the 
maps by Snell (|1981[) and Swade et al. (|1989a|). The ^^CO 



al. |1998| ). The mass within 10' of position (0,0) is 25 Mq for 
d =100pc. Alternatively, the mass contained within the con- 
tour r(200 /im)=6-10^'* is 28 A/q. The main uncertainty is in 



emission increases towards the south while C'^Oand ^■^CO(l- 
0) delineate the dense core close to the given centre position. 
The ^■^CO(l-O) emission extends towards the southwest, how- 
ever, where the ^^CO(2-l) emission peaks well outside the 
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Fig. 7. Maps of the integrated line areas: i2cO(l-0), C^*(l-0), 
i3CO(l-0), "C0(2-l), H13CO+(1-0) and DCO+(2-l). The 
dashed lines in each frame indicate the location of the DCO+ 
core 



core seen in C^^O. The distributions of Hi3CO+(l-0) and 
DCO+(2-l) are clearly different from the CO species. The 
distribution of H^'^CO+(1-0) is similar to what was seen by 



Guelin et al. (1982). However, due to the better resolution of the 



present study the emission area is better resolved and the map is 
more structured. There is a second emission peak south of the 
centre at position (1,-4). The morphology of the DCOh-(2-1) 
emission follows closely that of the H^'^COh-(I-O). The emis- 
sion is concentrated along a narrow, but clearly resolved ridge 
with some extension west of the centre position. The western 
extension is seen in both lines, although it was not visible in 
the Guelin et al. H^^CO map. In the centre, the DCOh- emis- 
sion seems to be more concentrated around the peak position. 
However, this is mainly due to the smaller beam size of the 
DCOh- observations. 

Maps of Ci*O(l-0), 13CO(1-0) and i3CO(2-l) emission in 
different velocity intervals are shown Figs. |^-|ic|. In the low- 
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Fig. 8. C^*O(l-0) integrated antenna temperature in four ve- 
locity intervals. The velocity intervals are given in the upper 
part of the panels. The range of contour levels and the step 
between consecutive contours are written at the bottom of the 
panels. 

est radial velocity bin, C^^O emission is seen both northeast of 
the position (0,0) and towards the western edge of the map. In 
the interval 2.0-2.5 km s^^, strong emission still exists in the 
western part of the map. The main core is visible in all three 
higher velocity intervals. Below 3 km s~^ the emission extends 
between the positions of the 100 fim and 200 /im maxima (see 
Fig. ||) while at the highest radial velocities the emission is con- 
centrated around position (0,0) only. The same general mor- 
phology and dependence on the radial velocity is repeated in 
the ^■^CO maps. The relative intensity of the ^■^CO lines is, 
however, higher in the south. As the radial velocity increases 
the emission maximum moves over the C^^O core from north- 
east to southwest. 

The DCOh-(2-1) and U^^CO+(2-l) in Figs. show 
similar behaviour At low velocity, the emission is concentrated 
in the south around position (2,-3) while at other velocities 
the maximum is reached close to position (0,0). There are also 
some differences, and H^'^COh- is relatively stronger around the 
emission peak in the south. 

4. Models of radiative transfer and column densities of 
molecular gas 

LTE column density estimates were computed based on the 
^^CO(l-O) and C^^O(l-O) observations. The derived excita- 
tion temperatures are between 7 and 9 K. Compared with the 
rest of the cloud, the excitation temperature seems to be sys- 
tematically 1 K lower in the northern part close to the column 
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Fig. 9. ^^CO(l-O) integrated antenna temperature in four ve- 
locity intervals. The velocity intervals and the contour values 
are indicated in the panels. 
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Fig. 11. DCOh-(2-1) integrated antenna temperature. The ve- 
locity intervals and the contour values are indicated in the pan- 
els. 



density maximum. However, the difference is small and we de- 
rive column densities assuming a constant excitation tempera- 
ture of Tcx=8.4K. The LTE column densities were calculated 
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Fig. 12. H^'^COh-(I-O) integrated antenna temperature. The ve- 
locity intervals and the contour values are indicated in the pan- 
els. 
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Fig. 10. ^'^CO(2-l) integrated antenna temperature in four ve- 
locity intervals. The velocity intervals and the contour values 
are indicated in the panels. 



according to the formula given by Harjunpaa et al. (1996) and 
are shown in Fig |l3[ In the conversion from C^^O to H2 column 
densities, a fractional C^^O abundance of 1-10^^ was assumed. 
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Fig. 13. LTE hydrogen column density as derived from 
the ^^CO(l-O) and C^^O(l-O) observations. The contours 
start at value N(H2)=5-10^^cm^^ and are drawn at inter- 
vals of S-lO'^^cm"^ (the thick contour line denotes the value 
N(H2)=M022cm^2-)^ 



corresponding to the quiescent cloud positions in Harjunpaa et 
al. ( |1996| ) and Frerking et al. ( |1982| ). 

From the derived column densities we obtain a value of 
~47Mq for the cloud mass within 10' of the centre posi- 
tion, or some 35 Mq for the mass contained within the contour 



A^(H2)=1.0-10^^cm-^ 

For comparison, the C^^O and ^^CO emission were mod- 
elled with spherically-symmetric model clouds where the ra- 
diative transfer prob lem w as solved with Monte Carlo meth- 
ods (see Juvela et al. 1997 ). While spherical symmetry presents 
an extreme simplification of true density distribution, the mod- 
els will still be more realistic than the assu mption s underlying 
the LTE calculations (see e.g. Padoan et al. 2000 ). The model 
clouds are assumed to be isothermal with radial density depen- 
dence of n ^ r~^-^ and with density contrast 20 between cloud 
centre and outer boundary. The free parameters of the models 
were the cloud size, central density and the kinetic temperature. 
The Ci*O(l-0), "CO(I-O) and "CO(2-l) spectra within 10' 
radius of the position (0,0) were compared with spectra calcu- 
lated from the radiative transfer model and the free parameters 
were adjusted in order to obtain the best fit. The gas distribu- 
tion is asymmetric with respect to the centre position and the 
^■^CO emission is displaced with respect to C^^O. Therefore, 
the results can not be expected to be very accurate. However, 
although the listed free parameters are not well constrained, the 
column density estimates should be more reliable. 

Assuming a constant kinetic temperature of lOK, the best 
fit was obtained with a model having a mass of 40 A/q . The 
mass estimate is not very sensitive to the assumed kinetic tem- 
perature. A similar model with Tkm increasing linearly from 
8.0 K in the centre to 15 K on the cloud surface resulted in a 



mass estimate of 34 Mq. The quality of the fit was, however, 
marginally better in the isothermal model. 

The spectral lines do not show signs of cloud collapse. 
Several ^^CO spectra have double peaked, asymmetric pro- 
files but comparison with optically thinner C^^O indicates 
that these are caused by separate emission components rather 
than self-absorption. Gaussian fits give average linewidths of 
-0.7kms-i for both i3cO(2-l) and Ci^O(l-O). The virial 
mass was estimated according to equation 

M^^ = i^fR, (1) 



where ctsd is the three-dimensional velocity dispersion, R 
the radius of the cloud and C is a constant depend- 
ing on the assumed model of the cloud. The value of C 
is 0.0508 km s^^ y^pc/MQ for a homogeneous sphere and 
0.0463km s~^ yJ pcjM^ for density distribution n ~ (see 
Liljestrom 1991 ). Assuming gas temperature 13 K and R ^ 10' 



the observed velocity dispersion corresponds to virial masses 
of 54 Mq and 44 Mq for the two models. In the presence of 
noise the FWHM derived from gaussian fits may, however, be 
biased and the average FWHM read directly from the spec- 
tra is smaller, ~0.5kms^^, lowering the virial mass estimates 
to 36 A/q and 30 A/q, correspondingly. The previous analysis 
ignores both the external pressure and the presence of veloc- 
ity gradients. The cloud is, however, approximately in virial 
equilibrium. Even though the cloud as a whole may not be col- 
lapsing it can still contain collapsing fragments. In fact, there 
is evidence of several pre-protostellar cores in L183 (Lehtinen 
et al. [2000| ). 

We will return to the modelling of the molecular line data 
in a future paper (Juvela et al., in preparation) where radiative 
transfer models will also be constructed for the far-infrared dust 
emission. 



4.1. Comparison of FIR and molecular line maps and column 
densities 

Fig. |l^ shows the distribution of FIR surface brightness rela- 
tive to the intensity of the C^^O(l-O) lines. The 200/xm inten- 
sity closely follows the C^^O emission. The peak positions do 
coincide and the second 200/im maximum in the southeast as 
well as fainter extensions towards south and towards southwest 
have corresponding features in the line map. The 200/im emis- 
sion extends, however, somewhat further in the north while the 
C^^O is relatively stronger in the southwest. With respect to the 
200 fim. peak the C^*O(l-0) maximum is displaced slightly to- 
wards the southwest. Further towards the southwest, the C^'^O 
is again well traced by the FIR emission. The second C^*0(1- 
0) peak at position (4,-2) coincides with the second 200 /^m 
peak. There is, however, a third 200 /im peak at (1,-6) which 
does not have a direct counterpart in the molecular line emis- 
sion. The peak corresponds to just one pixel above the diffuse 
background and could even be caused by a single extragalac- 
tic source. The lOO/^m peak traces the eastern extension of the 
C^^O distribution. Apart from this there is a clear lack of corre- 
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Fig. 14. The far-infrared maps of L183 at lOO^m and 
200/im with C^^O(l-O) contours overlaid. The contours are 
drawn between 0.4Kkms^^ and 1.6Kkms~^ at intervals of 
0.3Kkms-i 



lation with molecular line emission, especially within the C O 
core. 

In Fig. |15| Ci*O(l-0) and ^^€0(2-1) Une areas are com- 
pared with the 100/xm and 200/im surface brightness val- 
ues. Of the two far-infrared bands 200/im is clearly better 
correlated with molecular material. While on a logarithmic 
scale the relation between C^'*O(l-0) and 200 fj,m emission is 
approximately linear, the ^^CO(2-l) intensity becomes satu- 
rated beyond 5(200 //m)=50 MJy sr"^ Based on IRAS obser- 
vations Laureijs et al 1995 calculated the difference A/ioo = 
/(100/im) — /(60/im)/0.21. The 60 /im emission traces the 
diffuse outer parts of the cloud and the correlation with the 
^■^CO core was found to be better for A/ioo than for the original 
100 /im. The 100 /im observations are affected by this warmer 
emission and this causes the poor correlation. The colour tem- 
peratures south of the LI 83 core could also be biased due to 
diffuse material the presence of which is evident from the ^^CO 
maps (Fig. 

In Fig. 

Sect. 0) is compared with the dust optical depth at 200/im. 



16|, the C^^O(l-O) column density (see Fig. 13 and 



The figure confirms the good general correlation between the 
two. The C^^O column density distribution may, however, be 
slightly more extended in the SE. This was already seen in 
the C^^O(l-O) distribution. On the other hand, 200/im optical 
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Fig. 15. Ci8O(l-0) and i3cO(2-l) Hne areas as the func- 
tion of 100/im and 200/im surface brightness. In plots contain- 
ing the 200 /im surface brightness values the original data are 
smoothed to resolution of 100". In the case of 100 /im the cor- 
responding resolution is 70" 



depth is relatively larger immediately south of the (0,0) posi- 
tion. 

Nevertheless, the linear relation does hold only at lower 
column densities. A linear least squares fit was made taking 
into account the uncertainty in both variables. For points with 
T(200/im)<5-10-3 we get 

= (4.62 ± 0.16)10'' r(200/im) - (1.5 ± 1.0). (2) 

lO^'^ cm^^ 

Points at r(200/im)<510^^ were also discarded for this fit 
since the C^^O column density estimates corresponding to 
these points are no longer reliable. When T(200/im) exceeds 
'^510^'^ the C^'^O column density is significantly below the 
prediction based on lower optical depths. The C^^O emission 
could be reduced by depletion of gas phase molecules, a large 
drop in the excitation temperature towards the cold core or by 
the effect of optical depth. 



5. CO depletion in the cloud core 

The LTE analysis indicated C^^O(l-O) optical depths below 
one. In the Monte Carlo model, the optical depth towards the 
cloud centre averaged over the telescope beam was 0.56. These 
results were based on the ^^CO(l-O) and C^^O(l-O) observa- 
tions. Since the emission of the observed lines may come from 
different parts of the cloud, the determination of the optical 
depth remains uncertain to some extent. In the positions close 
to the cloud centre where C^*0(2-l) was observed LTE anal- 
ysis based on the two C^'^O lines gave somewhat higher esti- 
mates for the optical depth, t <1.0. The optical depth is insuf- 
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ficient to cause any significant saturation of the C^^O line and 
will not affect the accuracy of the column density estimates in 
Fig.[l|. 

Based on the LTE analysis, the excitation temperature was 
rather uniform over the cloud. Compared with other parts of 
the cloud the excitation temperature towards the cloud centre 
was lower by no more than 1 K. The column density estimates 
were based on an average value of 8.4 K. In Fig. |l^ the vertical 
arrow indicates the change in the estimated column density if 
excitation temperature is reduced from a value 8.8 K to a value 
of 7.0 K. It is clear that even a drop of ~2 K is insufficient to 
explain the relatively low C^^O column density estimates at the 
position of high values of T(200/xm). 

The horizontal arrow in Fig. |l^ indicates the shift in the 
value of T(200/xm) when a value of optical depth is calculated 
with rdust=12.5K instead of 12.0 K. Incorrect dust tempera- 
tures are, however, not likely to explain the relative drop of 
the C^^O intensity. Firstly, the error should be '^l K or more 
at 12K to bring the deviating points onto the fitted line even 
if the temperature were not changed outside the core. Errors 
in temperatures could be caused either by calibration or the 
background subtraction. If all 100/im surface brightness val- 
ues are multiplied by a constant, all T(200/Ltm) estimates would 
be scaled with a factor that has very weak temperature depen- 
dence. An incorrect background subtraction, on the other hand, 
mostly affects the high temperatures which correspond to low 
surface brightness values. As an example, let us consider two 
dust temperatures, 12 K and 15 K. After the background sub- 
traction, in L183 these typically correspond to 100 /im surface 
brightness values of 7MJy sr^^ and 3MJy sr^^. If we subtract 
O.SMJy sr^^ from the 100 /im measurements the colour tem- 
peratures would drop to 11.85K and 14.45 K and the corre- 
sponding 200 /im optical depth estimates would increase by 8% 
and 20%, respectively. The net change is only some 12% i.e. 
much less than the effect seen in Fig. p^ . 

Therefore, it seems that C^*0 depletion remains the most 
probable cause for the decreasing iV(C^^O)/T(200/im) ratio. 
In cloud cores the CO molecules freeze out onto the surfaces of 
the cold dust grains. Direct evidence of the process is provided 
by the ob servations of the molecular ice features (e.g. Tielens et 
al. 1991 ). As the result the properties of the dust grains are al- 
tered and this may explain some of the colour temperature vari- 
ations observed in LI 83. In the centre of LI 83 the visual extinc- 
tion exceeds 10™ and CO depletion has been detected in other 
clouds at similar extinction values. Kramer et al. (1999) and 



2.S 10' 



Bergin et al. (2001) have detected CO depletion in the cloud 
ICS 146. In that cloud the relative C^^O intensity decreases al- 
ready below Av=10™. Similarly in the cloud L977 a drop in 
C^^O intensity attributed to depletion takes place at visual ex- 
tinction Ay - 10™ (Alves et al. |1999| ). In L183 our models 
predicted central density close to 10^ cm^-^ for the spherical 



model with a radial density profile 



-1.5 



This density is 



probably exceeded in the cold core and e.g. in L1544 CO de- 
pleti on wa s associated with densities above cm^^ (Caselli 
et al. |1998|). 
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Fig. 16. a) Correlation between the dust optical depth at 
200 /im and the C^*0 column density derived from C^*0(1- 
0) observations, assuming a constant excitation temperature of 
8.4 K. The arrows indicate the shift if iV(C^*0) were calcu- 
lated with a ^2K lower excitation temperature, or if the optical 
depth were calculated using a 0.5 K higher dust temperature. 
The latter applies to the highest r(200/im) values correspond- 
ing to Tdust ^12K. b) The distribution of the 200 /im optical 
depth over the area of highest optical depth. The dotted contour 
indicates the level T(200/im)=210^'^. In (a) the points inside 
this contour are plotted with larger symbols. The solid contours 
show the distribution of the C^^O column density. The highest 
contour is at 1.6- 10^^ cm^"^ and the others are each a factor of 
two below the previous one. 



Based on Fig. |l6[ the depletion factor is '^1.5 in the centre 
of LI 83 where the dust temperature is close to 12 K. This is 
in perfect agreement with the findings of Kramer et al. ( 1999| ) 
who have derived the depletion factor for ICS 146 as a function 
of dust temperature. 
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6. DCO+ in the cold cloud core 

H^'^COh- and DCO+ are concentrated close to the 200/im max- 
imum. The peaks of the intensity distributions of both species 
are very close to the position of the continuum source, within 
~ 2'. While the 200/im emission continues to be strong to- 
wards the 100/im peak both H^^COh- and DCOh- are restricted 
to a narrow ridge running in the north-south direction. The 
distribution is therefore clearly different from either C^^O or 
200/im, although close to position (0,0) the 200/tm emission is 
also elongated in the north-south direction. More importantly, 
no H^^COh- or DCOh- emission is seen close to the 100/im 
maximum, while the area of strong C^*0 emission covers also 
that position. At the position (4.5', -2.5') we get 2a upper limits 
of O.IOK for H13COh-(1-0) and 0.20K for DCOh-(2-1) main 
beam temperature. 

Fig. |l^ shows the DCOh- distribution in relation to C^^O. 
Both peak in the same region around position (0,0) and the 
small extension from this position towards the west can be in- 
terpreted to correspond to a similar feature seen in C^^O. How- 
ever, at the exact position of the DCO+ maximum, the C^^O 
map shows a very slight depression, and in the south the two 
species are clearly anticorrelated, with stronger C^^O emission 
seen on each side of the DCO+ ridge. 

The 100/im emission peaks some 5' east of the DCO+ 
maximum and even in details the distribution is completely un- 
related to the DCO+ emission. The correlation between DCO+ 
and 200 [im emission is not very good either. Both peak close to 
position (0,0) but compared with the FIR emission the DCO+ 
maximum is shifted one arcminute to the west. Apart from a 
similar shift the correlation is good north of the centre position. 
In the south the distributions of DCO+ and 200 /im are differ- 
ent. As already mentioned, no DCO+ emission was seen close 
to the eastern FIR peak. The DCO+ map does not quite extend 
to the southern 200 /im peak at position ( 1 ,-6) but between the 
northern and southern 200 /tm peaks the DCO+ ridge goes over 
a region with relatively low FIR emission. 

In Fig. we compare the DCO+ distribution with the 
200/im optical depth. Due to the dust temperature gradient, 
r(200 /im) is shifted relative to the 200 /im emission and co- 
incides with the DCO+ distribution. Around position (0,0) and 
in the north, the correlation between r(200 /tm) and DCO+ is 
very good. In the south, the DCO+ follows roughly the dust 
optical depth distribution but the correlation is not perfect. The 
r(200 /tm) distribution is more extended SE of the centre posi- 
tion and the southern DCO+ peak at (1.5, -4) lies just north of 
the third T(200/im) peak, i.e. it does not follow dust distribu- 
tion. The southern part of the DCO+ ridge is outside the coldest 
dust core (Fig. 18 3). 

As the distributions of the DCO+ and H^^COh- are very 
similar, the same conclusions apply to the comparison of the 
H^^COh- distribution with CO species and FIR emission. At 
the position (0,0), the observed main beam temperatures are 
1.5 K and 0.79 K for DCO+(2-l) and Hi3cO+(l-0) respec- 
tively. DCO+ seems to be more concentrated in this core, but 
the difference is not significant when the difference in the beam 
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Fig. 17. Contours of DCO+(2-l) integrated line area on the 
map of integrated C^*O(l-0) line area. The DCO+ contours 
are at 0.2, 0.4, and 0.6Kkms~^. The dotted contours indicate 
C^^O(l-O) intensity levels from 0.5 to 2.0Kkms"i at inter- 
vals of 0.5 Kkms^^. The maps are presented at their original 
resolution. The beam sizes (46" for C^^O and 35" for DCO+) 
are indicated with circles at the lower left corner. 
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the 



size is taken into account. In the upper panel of Fig. 
DCO+ and H^'^COh- intensities are compared along the line go- 
ing through the two cores. The DCO+ values were convolved 
to the resolution of the H^'^COh- observations. In the main core, 
the H^^COh- is displaced to the north with respect to the DCO+. 
In the southern core, the relative intensity of H^'^COh- is higher. 

According to chemical models (Millar et al. 1989 ; Roberts 
& Millar |2000|) the ratio between DCO+ and H"CO+ can be 



used as a tracer of gas temperature. The reaction 



HD ^ H2D+ + H2 + AE 



(3) 



is strongly temperature dependent. At low temperatures the 
production of H2D^ is enhanced and it reacts with CO to pro- 
duce more DCO+. The abundance ratio also depends, however, 
on other factors such as cloud age, elect ron de nsity and the total 
abundance of neutrals (Anderson et al. 1999 ; Roberts & Millar 
2000|). This prevents the use of the ratio [DCO+]/[Hi3cO+] for 



direct temperature determination. 

The observed change in the intensity of the DCO+ and 
H^^COh- lines indicates a higher gas temperature in the south- 
ern core. This is consistent with the difference in dust temper- 
atures. Using the main beam antenna temperatures, the ratio of 
DCO+(2-l) and H"CO+(1-0) line ai-eas increases from ~1.0 
in the south to ^1.7 in the main core. The excitation states 
of DCO+ and H^'^COh- should be very similar. The isother- 
mal Monte Carlo model used to model the ^^CO and C^'^O 
lines was used to predict DCO+ lines assuming an abundance 
of 510"^^. Compared with the model in Sect. ^ the density 
needed to be increased almost a factor of two to produce lines 
Ta(J =1 — 0) ~0.8K, in which case the Une ratio was 
T(l — 0)/T(2 — 1) ~ 2. Assuming column density is directly 
proportional to observed intensity and using a value of 66 for 
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Fig. 18. a) DCO+(2-l) contours oti the map of 200 opti- 
cal depth. The solid contours show the DCOh- distribution and 
are drawn at 0.2, 0.4, and 0.6Kkms^^. The dotted contours 
indicate values MO"^, 2-10"^ and 4-10^^ of 200 optical 
depth, b) DCOh-(2-1) contours on the map of dust colour tem- 
perature. The solid DCOh- contours are drawn at 0.3, 0.5, 0.7 
and 0.9 Kkm^^. The dotted contours correspond to dust tem- 
peratures from 12 to 14 K at steps of 0.5 K. The maps are pre- 
sented at their original resolution. 



the ratio ^^C/^^C (Langer & Penzias 1993), we obtain values 
of 0.051 and 0.030 for the DCOh- to HCO+ abundance ratio in 
the two cores. The ratios indicate gas temperatures ^lOK or 
le ss, dep ending on the set of reac tion constants used (M illar e t 
al.|l989[ Roberts & Miliar- |2000l see also Loren et al. |l990|). 



The difference between the northern and southern cores corre- 
sponds to a temperature change of at least 2 K. 



7. Comparison of optical extinction with FIR dust and CO 
column densities 

In Fig. ^ we plot Ay{ as a function of as derived from the 



starcounts (see Sect. 2.3). The dots denote individual extinc- 



tion values obtained by averaging the counts in each grid point 
over a circle with radius 1.7'. In the averaging the weights were 
proportional to the area inside the circle which was, however. 
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Fig. 19. a) The DCO+(2-l) (solid Hne) and Hi3cO+(l-0) 
(dashed line) intensities along the line going through the two 
cores detected in these lines, b) DCO+(2-l) contours (solid 
lines) over-plotted on the Hi^CO+(l-0) map. The DCO+ con- 
tours levels are the same as in the previous figures. The dotted 
lines indicate H"CO+(1-0) values 0.2, 0.3 and 0.4Kkms-^ 
The DCO+ data are convolved to the resolution of H^^COh- 
observations. 



only slightly larger than original 2.8' x 2.8' square used in the 
star counts. 

We first calculated the average number of _R-band stars in 
those 2.8'x2.8' squares with a given number of _B-band stars. 
The average counts in these areas were transformed into ex- 
tinction values. The relation with statistical errors is shown in 
Fig. ^ (dashed line). Linear least squares fit gives a relation 

= (0.42 ± 0.02) Ab- The fit takes into account uncertainty 
in both Ab and Ar. If the intercept were included in the fit 
its value would be only ^ 0.03. If we calculate the average of 
the values at all grid points with a given Ab (instead of 
first calculating the average number of stars) the derived rela- 
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Fig. 20. The relation between B-band and i?-band extinction 
as derived from the starcounts. The lower points with error- 
bars are obtained by first averaging the i?-band counts in areas 
with given number of i?-band stars. The dashed line is the least 
squares line fitted to these. The upper points denote averages 
of An for given Ab and the dots correspond to star count av- 
eraged over circle with radius 1.7'. The solid line is the least 
squares fit to these points and has a slope of 0.48 



tion becomes steeper, with a slope of 0.48±0.02. The value is 
still somewhat lower than the value ^^0.56 for the extinction 
law with i?v — Ab/ Eb-y=3.1. The fit is shown as a solid 
line in the figure. The least squares line fitted directly to (Ab, 
A^) points gives also a slope of 0.48. Considering the possi- 
ble systematic errors of ^10% of our extinction values, due to 
uncertainties in the calibration of logN{m) vs. m slopes, our 
A^/Ab ratio is compatible with the standard reddening law 
(i?v — 3.1). On the other hand, it differs significantly from 
the value of ^r/^b=0.66 valid for the case of i?v = 5.0 as 
derived for the "outer-cloud dust" 



19901 



' (Mathis et al. 

In the centre of L183 there were a few starcount squares 
where no stars were seen and for Fig. ^ the extinction was 
calculated using a value of 0.5 stars per square. The values, 
A^ «5™, still only represent a lower limit of the true extinc- 
tion. In the most opaque area of 95 square arcminutes, four 
stars were detected in the i?-band and this translates into an ex- 
tinction of 5.3'". In the _B-band, no stars were seen in this area 
while a single star would correspond to an extinction of 6.8™. 
Therefore, it is likely that at the very centre the extinction is at 
least Ay -10"\ 

Next j4r is compared with the hydrogen column density es- 
timate derived from the C^®O(l-0) observations (Fig. |l3|) and 
the 200 iim optical depth (Fig. The data are convolved to 
the resolution 2.8' corresponding to the grid size used in the 
starcounts. The results are shown in Figs. ^ and |2[ 
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Fig. 21. a) i?-band extinction vs. C^^O column density, b) R- 
band extinction contours drawn on the C^^O column density 
map. The contours of extinction are drawn between 1™ and 
5'" at intervals of 1™. The C^^O column densities have been 
convolved to the resolution of the extinction map, '^2.8'. 

The C^^O column density and the extinction are clearly 
correlated. The C^^O depletion seen in Fig. |l^ is not visible, 
however, as the centre region where no i?-band stars were de- 
tected is not included in the figure. This limits the dynamical 
range available. Furthermore, at low column densities the noise 
increases rapidly due to the decreasing signal-to-noise ratio of 
the C^^O spectra. The maximum extinction corresponding to 
one star per grid pixel area was ~4"\ For the plot the star 
counts were again averaged over a circle with radius 1 .7' but 
points with A^ >4™ are omitted in the figure. The least squares 
hne is 

An/mag = (3.4 ± 0.2)1Q-^'^ N{C^^O) + (0.1 ± 0.1). (4) 
With the earlier relation between ^Ir and Ab and the assumed 
value of i?=3.1 this translates to 7V(CiSO)/cm-2=l.9.l0i4Av- 
4- 10^'^. The slope is similar to the value found in regions 
of high extinction in Taurus and Ophiuchus by Frerking et 



al. (1982) but steeper than the corresponding value for the 
< Ay < 4™ where the value in Taurus 



lower range 2" 

was 0.7-10^'* cm^^ mag^^. The slopes determined in L977 and 
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Fig. 22. a) i?-band extinction vs. 200 fim optical depth, b) R- 
band extinction contours on the map of T(200/im) which has 
been convolved to a resolution of 2.8'. The extinction contours 
are drawn between 1'" and 5™ at intervals of 1"\ 



ICS 146 for Ay ^10'" (Alves et al. 1999) are rather similar. 



'^2-10 cm^ mag^ . Also Harjunpaa et al. (1996) found a 
similar slope in Corona Australis while in the Coalsack the 
value was only ~1.2-10^"'cm^^ mag^^. Note that if we assume 



standard extinction (Rieke & Lebofsky 1985) to transform 
directly to Ay our slope increases to 2.2-10^'* cm~^ mag~^. 

Within the statistical uncertainty of the starcounts, the ex- 
tinction map agrees well with r(200/xm). The dots in Fig. ^ 
correspond to the star counts on the original grid. The filled 
squares with error bars are averages for each discrete value of 
Ar. The least squares line based on the interval 1.0'" < A^ < 
4.0" is 

An/m&g = (2.2 ± 0.3)10V(200/xm) + (0.2 ± 0.2). (5) 
The relation can be used to derive an extrapolated extinction 
value for the centre region where no stars were detected. In 
the centre there are several 200 /^m pixels with r(200/im) ~ 
6 • lO^'^. According to the previous relation, this corresponds 
to extinction « 13™ and visual extinction of the order of 
17™. This clearly exceeds the previous estimates, and is based 



on the sharp increase of the 200 /im optical depth in the cold 
cloud core. 

Assuming a standard extinction law (Ry=3.1) between 
V- and _R-bands, from Eq. |] we obtain T{200fira)/ Ay = 
3.4 • lO^^'mag^^. For the cloud centre, the ratio increases to 
3.8-10^'' mag^^ but is still below the value of 5.3-10~''mag^^ 



found in the Thumbprint Nebula (Lehtinen et al. 1998). In 
L183, the comparison between 100 /xm optical depth and vi- 
sual extinction gives a ratio 2.1-10~^mag^^ which is similar to 
the value in the Thumbprint Nebula. 

For a value of r(200/im) = 1.0 • 10^"^, the pre- 
vious relations between T(200/im) and Ar and between 
Ar and C^*0 column density predict a column density of 
N(C^^O)=6.610^^cm^^. Directly using the relation between 
r (200/im) and N(C^^O) (Eq. ||) we would obtain a lower value, 
4.5-10^^ cm~^. One must remember, however, that Eq. H was 
determined over a much larger range of optical depths and con- 
sequently that relation may akeady be affected by the C^^O 
depletion. 



8. Conclusions 

The main conclusions drawn from the comparison of far- 
infrared data and molecular Une observations in L183 are: 

- The 100/im and 200/im emission clearly have different 
distributions. 

- The 200/im maximum indicates the position of a core 
where either the dust temperature becomes very low or the 
dust properties have changed. The colour temperature de- 
termined from the 100/im and 200/im observations with 
an emissivity law of is ^12 K. The dust optical depth 
peaks at the same position. 

- Of the molecular lines, the optically thin C^'^O best traces 
the dust distribution. Although better correlated with the 
200 /im surface brightness it shows strong emission at the 
locations of both the 100 fim and 200 maxima. The op- 
tically thick ^^CO and ^■^CO lines peak south of the dust 
cores. 

- C^*0 is depleted in the cloud centre. The depletion factor 
is ^^1.5 in the core where dust temperature is close to 12 K. 

- Relative to each other, the DCOh- and H^^COh- lines show 
very similar distributions. Their distribution is completely 
different from any of the observed CO lines but a good cor- 
relation does exist with the 200 /im optical depth. 

- The estimated mass within 10' of the (0,0) position was 
25 Mq, based on the 200 /im optical depths and ^40 Mq 
based on the C^^O column densities. 

- The 100/im surface brightness is not well correlated with 
the extinction. A linear relation does exist between Ar and 
dust optical depth at 200 /im. Based on this correlation we 
derive for the cloud centre an extinction Ay ^ 17™. 
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